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The definition, properties, and application of an extended correlation 
function of two random variables involving two common parameters are 
described and applied to mobile radio systems. The correlation functions of 
a predetection diversity combined signal (using a scheme of phase equalizing 
by multiple heterodyning) and of a directional antenna array signal are 
derived with the help of the extended correlation function. 

These correlation functions can be used to determine parameter values 
giving minimum correlation between two signals desirable for diversity 
systems. One can also obtain the power spectra by taking the Fourier trans- 
form of these correlation functions. Thus extended correlation functions 
promise to be useful. 

I. INTRODUCTION 

If two random variables depend on only one common parameter, 
such as time or distance, conventional correlation formula can be ap- 
plied to the two variables. However, if both of these variables involve 
not one but two common variable parameters, then the correlation 
formula found in the current literature is limited.* Since such cases occur 
in some of our mobile radio problems, as we discuss later, we need to 
define an extended correlation function and outline its properties and 
applications. 

II. DERIVATION OF AN EXTENDED CORRELATION FUNCTION OF TWO 
RANDOM VARIABLES INVOLVING TWO COMMON PARAMETERS 

A conventional normalized correlation function of two random 
variables r t and r 2 , both of which are functions of one parameter 

* Prior to acceptance of this paper for publication, the author was advised that 
a similar concept was discovered independently by A. Papoulis in his recently 
published book. 1 
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d[that is, r,(d,) and r 2 (d 2 )] can be expressed as' 

/?i 2 (^i i d 2 ) — W1W2 



.2 



p,M • A) - 



0-10-2 



where m's are the mean values, <r 2 's are the covariances, p, 2 (d, , da) is 
in the range £ | p 12 (d, , d 2 ) | ^ 1, and Bu(dj , d 2 ) = (r 1 (d 1 )r 2 (d a )>B» 
is the correlation function. + 

Supposing a random variable r x (Di ; d,) is a function of two parameters 
D, and d, , and another variable r 2 (D 2 ; d 2 ) is a function of two param- 
eters D 2 and d 2 ; the normalized correlation functions of these two 
variables can be deduced from equation (1): 

p, a (Z>, , D 2 ; rf, , d 2 ) 

fii 2 (£), , Z) 2 ; rfi , rf 2 ) ~ ''M"i 

fr,(D, , d,)r,(Z) 2 , d,)). v - (r,(i), , rf,)Ur,(Z), , rf 2 )) nv 

- [(^(d, , *)>.. - p; , d,)>a*[«(A ; rf 2 )>..v - <r 2 (D 2 ; iwsj 1 

(2) 

If the problem we are dealing with is a stationary random process for 
both of the parameters D and d, then 

72,2(2), , D 2 ;d, , d 2 ) - i2 I2 (Z>, - D 2 ;d l - d 2 ), 

(r k (D h , d k ))„ = (r,(0, 0)>„ - m k , 

(rl(D k , d k ))„ - ml = (rl(0, 0)> uv - ml = <r\ , 

where k = 1,2. Now ra* and a k are constants and we may let D = Z), — 
Z) 2 and d = di — d 2 . Then equation (2) becomes 

, n R l2 (D;d) - m 1 m 2 ,. 

o x a 2 

We call p, 2 (D; d) a normalized extended correlation function of the first 
kind. Also we note that p, 2 (Z); d) in equation (3) is always smaller than 
p, 2 (0; 0) which is equal to one: 

Pi.0>;d) S p,.(0;0) = 1. 



t The terras "correlation function Ri*(cL, d 2 ) and normalized correlation func- 
tion pia(di, da)" are adopted from Ref. 3, p. 59. 
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When the difference D is equal to zero, then 

p 12 (0; d) = p l2 (d). 

Now we should illustrate and extend equation (3). We are going to 
find an extended correlation function of the first kind from a function 
e(D; di , d 2 , d 3 , • • • , d m ) where all the d's are function of D and another 
parameter a [that is, d t (D, a) for i = 1, m], then 

R,(D; dt , da , d a , • • • , d m ) 

m <e[0;d 1 (0,0) > d a (0,0),d,(0 l 0), •••) 

■e(D; d^D, a), d 2 (D, a), • ■ • , d m (D, a)]) av , 

and the normalized correlation function can be derived from equation 
(2) as 

, n . , ,x R t (D\d x , • ■ • , d m ) - in, 

0". 

where 

m ( = <e[0; ^(0, 0), rf 2 (0, 0), d a (0, 0), • • • , d m (0, 0)]>, 1V 

a\ = (e 2 [0; d,(0, 0), d 2 (0, 0), rf 3 (0, 0), • • • , d m (0, 0)]> av - m\ . 

If we consider the case d { {D, a) is a constant for all D and a itself is 
a constant, then we may assign a new symbol R t {D \di , d 2 , • • • , d m ) 
which can be expressed as 

R t (D \d l ,d 2 , ■■■ ,d m ) 

= <e(0; d x ,d 2 , • • ■ , d m )e(D; d t , d 2 , d 3 , • ■ • , d n )) av . 

R,(D \ d t , d 2 , d 3 , ■ ■ • , d m ) is a correlation function under a condition 
that all dj , d 2 , d 3 , • • • , d m are constants. The normalized correlation is 

Pl (D | d> , r/„ ,.-.)« R.V>\d*,d a ,y,d m )-m\ (4) 

where m t and o-« have been defined previously. We call equation (4) 
the normalized correlation function of the second kind. As we will show 
in the Section III, the extended correlation function of first kind p i2 {D; d) 
and the extended correlation function of second kind p i2 (D \ d) can be 
used to obtain the correlation of signals from two diversity scheme re- 
ceivers easily. 

In order to give physical meaning to these functions, let us consider 
the following two cases. Suppose that two base-station multibranch 
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diversity receiver arrays are separated by a distance D. The antenna 
spacing between branch-antenna elements for the first array is d, and 
for the second array is d 2 . Both receivers simultaneously receive the 
signal from a distant mobile radio unit. We would like to determine 
the values of d v , d 2 , and D to obtain the least cross-correlation de- 
sirable for the best diversity reception of these two received signals. 
The extended correlation of first kind p l2 (D; d x , d 2 ) may be used in 
this case. 

The second case assumes that a mobile radio multi-branch diversity 
receiver array, with given uniform antenna element spacing d, moves 
along the street with a constant speed V. The autocorrelation of a 
signal e, received by the mobile receiver, can be obtained from the ex- 
tended correlation function of the second kind p,(D I d). Alternatively, 
we can also consider a multielement directive antenna instead of the 
diversity scheme. In this paper, we only treat the latter case. The 
former case can be solved following the same technique. 

III. APPLICATION TO MOBILE RADIO PROBLEMS 

3.1 Derivation of the Correlation Function of a Signal Received from a 
Predetection Diversity Combining Receiver 

A multichannel predetection diversity combining system is a scheme 
for bringing a number of RF carriers to a common phase by means of 
multiple heterodyning. Then a linear combiner at the IF frequency is 
used to sum the individual channels. 4 ,B A signal received from this system 
is called a predetection diversity combined signal. 

Suppose that a signal consisting of multipath vertically polarized 
waves is received by an M-branch predetection combining mobile re- 
ceiver with a il/-antenna space diversity array. The 7l/-antennas are 
spaced by d t , d 2 , ■ ■ ■ d M respectively from an arbitrary common point. 
After the array has moved a distance D, the received signal e, which is 
the sum of the M individual signal amplitudes received from M indi- 
vidual antennas, can be expressed as ' 

e(D; d, ,d 2 ,d 3 ,d 4 , • • • , d u ) = r,(D; d.) + r 2 (D; d 2 ) + r 3 (D; d 3 ) 

+ •■■ +r M (D)d at ) 

= !>„,(£>; rf,„), (5) 

i 

where all r m are functions of distance D and antenna spacing d m (see 
Appendix A). For a mobile radio signal, 7,8 or a long range fading signal, 
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the r m are usually Rayleigh distributed. Suppose that all d's are con- 
stants, then the autocorrelation function of the signal given in equation 
(5) as a function of the separation distance D is an extended autocorre- 
lation function of second kind which can be expressed as 

R t (D \d l ,d 2 , d 3 , di , • • •) 

= (6(0; d 1} d 2 , d, , d 4 , • • ■)«(£; d, ,d 2 ,d 3 , d 4 )). ¥ 



2>.(Z>;0 



Er.(0;W 

1 

- (E Er B (0;iWD;<y 

\m=l n=l 

= E E «-(£;«*. - d n ). 



(6) 



Using equation (3), this can also be written 

R ( (D \d, ■■■ d u ) = Pt (D | d 1 , ■ ■ ■ , d„)(0 + m\ 
where 

<x 2 = (e 2 (0;d, , ••• ,d M )) - m\, 



(7) 
(8) 



<6 2 (0;rf, , ••• ,d A/ ))av = ^2> m (0;dJ 

= E EW0;4K(0;a 
i i 

= E E^ mn (0;c/,„ - d B ). 



(9) 



Substituting equation (S) into equation (7), and combining equations 
(6) and (7), we obtain 

E iLR,„ n (D;d m - d„) - m\ 
P ,(D | d x , ■ ■ ■ , d M ) = ^-^ 



E ZRU0;d m - d n ) - m\ 

1 1 

The terms R mn (D; d m — d n ) can be found from equation (3); 
R mn (D; d m - d n ) = Pmn (D; d m - d n )a m a n + m m m n 
and 



(10) 



(n; 
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m\ = (e(0; d, ■■■ d m ))l = (l> m (0; d m )) 
= [£<r m (0;d m )) av ] 2 

At M 

= EE m m m n . (12) 

i i 

Hence the correlation function of equation (10) becomes, assuming 

M M 

E HpUD;d m -d n ) 
Pt (D I d t , ■ • • , d M ) = =^-^ (13) 

E EUQ;i-« 

m=l n=l 

If all spacings between two adjacent antennas are equal, then d m — 
d n = (m — n)^! where di is the distance between two adjacent antennas. 
We may let d = d x , and simplify the notation of equation (13) to 

M M 

E E pUD; d) 

pXD I d) = ^w (14) 

E Ep m n(0;d) 

m=l n=l 

Equation (14) shows that a normalized autocorrelation function of 
an M-branch predetection combined signal is a normalized extended 
autocorrelation function of second kind in terms of all individual nor- 
malized correlation functions between branches. We notice that 

Pt (D | d) ^ Pt (0 | d) = 1, (15) 

and as stated in Section II 

Pmn (0;d) = Pmn (d). (16) 

We may also realize that 

Pl2 (Z>; d) = P23 (D; d) = p 3i (D; d) = • ■ • 
and 

Pl3 (D; d) = p 24 (D; d) = Pa5 (D; d) = ■ • ■ . (17) 

Hence, equation (14) can be further simplified as 

pXD I d) 

_ i¥p 11 +(M-l)(p 12 +p 21 ) + (M-2)(p 13 +P3 l )+---+Pu/ + PAn 
" i¥p° 1 +(M-l)(p° 2 +p 21 )+(M-2)(p° 3 +P3,)+---+P?A/ + P°n ' 

(18) 
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where p mn «- p mn (D; d) and p°„(0; d) used in equation (18) are for sim- 
plicity (p mn and p^,„ are derived in Appendix A). If we let the antenna 
.spacing d/X = 0, then p t (D | 0) from equation (18) represents the 
correlation function of two single-branch signals 



Pl (D | 0) = JXfiD) 



(19) 



which agrees with that in Ref. 6. 

Several numerical calculations have been carried out for the following 
example: Two four-branch diversity receivers, each of them with fixed 
antenna spacing d/X = 0.5 or d/X = 1.0, are mounted on the roof of 
the mobile unit, as shown in Fig. 1. These two receivers are separated by 
a distance D/X (D/X varies from to 4) for two cases, a = 0° and a = 
90°. The calculations of the extended correlation function p ( (D \ d) of 
these two signals, obtained from their respective receivers when the mo- 
bile unit is moving, are shown in Figs. 2 and 3. Both figures indicate the 
values of D/X which give the least correlation between two signals. We 
also note that the correlations at a = 0° are higher than that at a = 90°. 
Figures 2 and 3 can also represent the auto correlation of a signal re- 
ceived from a single four-branch diversity receiver which has its antenna 
spacing d/X = 0.5 or 1.0 and moves on a street with a constant speed 
V(D = Vt). The power spectrum of such a signal can be obtained by 
taking the Fourier transform of its autocorrelation function. 

3.2 Derivation of the Correlation Function of a Signal Envelope Received 
from a Directional Antenna Array 

Signal reception from a directional antenna array with M antenna 
elements has been also suggested as a means of overcoming multipath 
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Fig. 1 — Coordinate system of a M -branch diversity mobile radio receiver 
(M = 4 branches). 
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D/X- DISTANCE WHICH A FOUR - BRANCH DIVERSITY RECEIVER MOVES 

Fig. 2 — Normalized autocorrelation function of a four-branch diversity re- 
ceiver moving at a = 0°. 

fading in mobile radio propagation. 10-12 The derivation of the correlation 
function of this signal envelope is as follows. 

Suppose that the same kind of signal which consists of multipath 
vertical polarized waves as mentioned in Section 3.1 is received by a 
directional M -antenna array. The M antennas are spaced by c/j , d 2 , 
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0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

D/X -DISTANCE WHICH A FOUR -BRANCH DIVERSITY RECEIVER MOVES 

Fig. 3 — Normalized autocorrelation function of a four-branch diversity re- 
ceiver moving at a = 90°. 
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d* ) • • • i d M respectively from an arbitrary common point. After the 
antenna array is moved by a distance D, (see Fig. 4) the received 
signal envelope e which is the amplitude of the sum of M individual 
signals can be expressed as 13 



e(Z);d! ,d 2 ,d 3 , •■■ , d M ) 
- \s l (D;d 1 ) + s 2 (Z);tf 2 ) + 



+ s u (D;d M ) 



= | X(D; d, ,d 2 , ■■■ , d„) + jY(D; d x , d 2 , ■ ■ • , d M ) |, (20) 

where s m is a complex variable which represents the amplitude and the 
phase of an individual signal. X and Y are the real and imaginary parts 
of the total signal. 

If the spacings between adjacent antennas are equal, then antenna in 
and antenna n are separate dyb d m — d n = (m — n)d. Therefore X and 
Y of equation (20) are functions of D and d only. Suppose that all d's 
are constants, the autocorrelation function of signal envelope e can be 
obtained by using the equation: 14 



P.(D | d) = 



(X,(0; d)X 2 (D; d))l + (Z,(0; d)Y 2 (D; d))l 
<X?(0;c/)>L 



(21) 
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Fig. 4 — Coordinate system of a broadside directional antenna array (M = 8 
elements). 
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provided X and Y are gaussian variables, all <X m ) av are zeros, and 
(X*) av and <F^> av are equal, where m = 1 or 2. These facts are shown 
in Appendix B. If the antenna spacing cZ/X = 0, then p t (D | 0) actually 
represents the correlation between two single-branch signals, which 
agrees with equation (19) and Ref. 6. 

The normalized correlation function of a signal received from a 
broadside directional antenna array is 



pXD \d) = \ 



where 



Z £ [JMJ + /.(BO + Jo(A 2 ) + J (B 2 )]\ 

m=l n=l J 



r k k 

Lm=l n-1 



[J (A ) + J (B )] 



(22) 



M 
K = — for Mis even 



2 

M + 1 



M for is odd, 



and A i , Bi , A 2 , and B 2 are shown in equation (48). A and B are 
shown in equation (49). 

Several numerical calculations have been carried out for the following 
example: Two eight-element broadside antenna arrays, each of them 
with fixed antenna spacing d/\ = 0.5 or d/\ = 1.0, are mounted on 
the roof of the mobile unit. These two arrays are separated by a distance 
D/\ (D/X varies from to 4) for two cases a = 0° and a = 90°. The 
calculations of the extended correlation function p ( (D \ d) between two 
signals received from their respective arrays when the mobile unit is 
moving are shown in Figs. 5 and 6. Both figures indicate the values of 
D/\ which have the least correlation between two signals. The extended 
correlation curve of d/X = 0.5 is quite different from that d/\ = 1.0 in 
both figures. The curve of d/\ = 0.5 in Fig. 5 shows that the high 
correlation and low correlation are about 0.25X apart; however, this 
phenomenon does not appear for d/\ = 0.5, but rather for d/\ = 1.0 
in Fig. 6. It can be explained as follows. For the directional antenna 
array with spacing d = X/2, most of the energy is contained in the two 
major broadside lobes, while for the directional antenna array with 
antenna spacing d = X, most of the energy is contained in the two major 
end-fire lobes. As the vehicle moves, strong standing waves may occur 
when the major antenna lobes he in line with the motion of the vehicle, 
such as for the case a = 0° and d = X/2; or the case a = 90° and d = X. 
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0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

D/X-DISTANCE WHICH AN EIGHT-ELEMENT DIRECTIONAL ANTENNA ARRAY MOVES 



Fig. 5 — Normalized autocorrelation function of an eight-element directional 
antenna array pointing at a = 0°. 

The autocorrelations obtained from these standing waves, then, be- 
come oscillatory in nature, as we would expect. 

Figures 5 and 6 can also represent the autocorrelation of a signal re- 
ceived from an eight-element broadside antenna array which has its 
antenna spacing d/X = 0.5 or 1.0 and moves on a street with a constant 
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0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

STANCE WHICH AN EIGHT- ELEMENT DIRECTIONAL ANTENNA ARRAY MOVES 



Fig. 6 — Normalized autocorrelation function of an eight-element directional 
antenna array pointing at a = 90°. 
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speed V (D = Vt). The power spectrum of such a signal can be obtained 
by talcing the Fourier transform of its autocorrelation function as we 
mentioned in Section 3.1. 

IV. CONCLUSION 

The derivation of a general correlation function of two random 
variables, each of them involving two parameters, has been obtained. 
The terms "extended correlation function of first kind" and "extended 
correlation function of second kind" have been defined. The application 
of the extended correlation function is demonstrated. The correlation 
function of a diversity signal and the correlation function of a directional 
antenna array signal are derived with the help of the extended correlation 
function in this paper. Several numerical calculations have also been 
carried out. From these correlation functions we can obtain the least 
correlations between two signals under certain circumstances. Also, 
we can obtain the power spectra by taking the Fourier transform of 
these correlation functions. Thus, it seems likely that these functions 
will find general application. 
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APPENDIX A 

Finding Normalized Cross Correlation Functions From Individual 
Branch Signals of a Predetection Diversity Combining Receiver 

It is easy to show that the signal from branch m in equation (5) is 

N 

r m = £ A » ex P [+## cos (0« - «) + Km - l)0d cos U ] 

i 

-IX. + JF.I, (23) 

where 

A u = R u + jS u , 

AT 

X m = £ #„ cos <f> u + S u sin U , (25) 

u-l 

N 

Y m = Y, S u cos 4> u - R u sin <f> u , (26) 
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-0„ = 0D cos (0 U - a) + (to - 1) &d cos 0„ (27) 

(# u and S u are independent gaussian amplitudes with zero mean and 
unit variance). The diversity receiver is located at the point (D, a) in 
polar coordinate. The distance D, the angle a, and the arrival of uth 
wave at angle 0„ are shown in Fig. 1. We assume the N waves are uni- 
formly distributed in angle. Now we can average the product of two 
components of two branches — branch m and branch n — as 

<X.(D,)Z.(Z) 1 + D)>.„ 

= (X m (f>)X n (D))„ 

= NE{ cos [BD cos (0„ - a) - (m - n)0d cos 0„]} 

= N[J Q (a)J (b) - 2J 2 (a)J 2 (b) + 2J 4 (a)J 4 (6) 

- 2J,,(a)J,Xb) + ■■■] 

= N.Jvia 2 + 6 2 ) 5 , (28) 

where 15 

a = @D cos a - (to - w)j8d, (29) 

6 = 07) sin a, (30) 

a 2 -|- b 2 = (BD) 2 + (to - n) a (0d) 2 - 2(to - n)0 2 Dd cos a, (31) 

and 

(X m (D l )Y n (D l + Z»> av 

- (X a (0)Y n (D))~ 

= NE \s'm [0D cos (6 U — a) — (m - n)0d cos U ]} 

- 0. (32) 
Also 

<X m 2 (Z>,)> = (J'^DO) = tf. (33) 

Substituting equations (28), (32), and (33) into the following equa- 
tion 6 ' 1 ' 

,_ « . (X,„(0; d)X„(D; d))l + (X,„(0; d)F.(J); d))U „,. 

Pmn (D, d) * vaoTd)Z (34) 

Then we obtain the final result 

p mn (D; d) = J 2 o [(0D) 2 + (to - n) 2 (f3D) 2 - 2(to - n) 2 Dd cos a] J 

(35) 
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and 

Pmn (d) = Pmn (0; d) = Jl[(m - n)0d}. (36) 

We also can show the relations 

Pl2 = p23 = P34 ! P21 = P3 2 = P43 » 

Pl3 = P24 — P35 i P31 = P42 = P83 > 

Pmn 5^ p„ m fOT W, 7± U, 

Pmn (D;d) = Pnm (-D;d). 

APPENDIX B 

Finding a Normalized Correlation Function From a Real Part and an 
Imaginary Part of a Signal Received From a Directional Antenna Array 

It is easy to show that a signal consisting of N multipath vertical 
polarized waves received from an equal-spaced directional antenna 
array at a distance D from a reference position is 

E,(D) d) = 2 4 U {1 + exp 0» + exp (j2$) 

u-l 

+ exp (jdrfi) + • • • + exp [j(M - 1)^]} 

•exp(jj8D cos 0„), (37) 

where A u was defined in equation (24), 

\p = fid sin (a — U ) + 8, 
d is antenna spacing between two antennas, 
M is the number of elements, 
a is the normal direction of the array, 
5 is the relative phase between antennas, 

D is the distance measured from the coordinate origin to the center 
position of antenna array. (The center position of the antenna 
array is assumed always on the axis, that is, at the position 
(D, 0).), and 
6 U is the angle of arrival of the uth wave and is assumed to be uni- 
formly distributed. 

The coordinate system of a directional antenna array is shown in 
Fig. 6. Since the spacings between antennas are equal, we can let the 
phase refer to the center point of the array. Then equation (37) can be 
simplified by combining the first term and the Mth term, the second 
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term and the (M — l)th term and so forth. 13 The result becomes 

N 

E.(D;d) = £ A u exp (j0D cos 6 U ) 

.[2cos(^^) + 2cos(^^) + ... + 2 Q ] ) 

(38) 
where 

Q = 1 if M = odd 

= cos (\\p) if M = even. 
Equation (38) can be separated into a real part and an imaginary part as 

E,(D; d) = X + jY 
and 

e(D; d) = | E t (D; d) \ = (A' 2 + Y 2 )\ (39) 

where 

K N 

Z = 2EE [fi» cos WD cos 0„) - S u sin (/3Z> cos U )] 

m- 1 « - 1 

/M + 1 - 2m 
• C0S l 2 

= 2 £ * m (40) 

m-l 

K N 

F = 2 £ £ [# u sin WD cos 6 U ) + £„ cos (/3Z) cos U )] 

'M + 1 - 2m 



cos 



m-l 




where 






rr M ., ,, . 

K = — if M is even 




M + 1 ., . , . , , 
= — - — if M is odd. 



(41) 



(42) 
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Since R u and S u are independent gaussian variables, it is easy to 
realize that all x m and y m are independent gaussian variables. Hence 
X and Y are also gaussian variables. The mean values of X and Y are 
zeros, and the mean squares of X and Y are the same. Therefore equa- 
tion (21) can be applied. The following term in equation (21) can be 
replaced by 



(X,(0; d)X,(D; tf)>„ = 4| E .r„,(0; d) E x,„(D; d) 
= 4 E Z(^(0;rf).x„(D;r/)> nv . 



(43) 



The term <Xi(0; d)F 2 (Z>; d)> also can be obtained, and is equal to equa- 
tion (43), by replacing x n by y n . Then equation (20) becomes 



P.(D/d) = 

[~E EWOiW^iC 



E E<«.(0;d)y.(D;«0>« 



_m=l n=l 



|~E EW0;W0;C 

Ltb-I n=l 



J- , (44) 



n-l n=l 

where K is shown in equation (42), and 
(x m (0; d)x n (D; d)) av = ^ (cos (0D cos 0„) • {cos [(Af + 1 - m - n)fl 

4- cos [(m — n)^]}) BV . 
(.c m (0; d)y n (D; d» BV = ^ (sin (/3D cos 0,,)- {cos [(M + 1 - m -n)f] 



and 



-f- cos [(m - n)^]}) av , 



^ = /3cZ sin (a - U ) + 5. 



(45 



Now we may consider only a broadside directional antenna array, 
that is, 5 = 0. Then the following terms can be derived: 15 

(cos (a cos 0„) • cos [b sin (a — d u )]) RV 

= £(cos [(a 4- b sin a) cos 6 U - b cos a sin U ] 

+ cos [(a — b sin a) cos U + 6 cos a sin 0„]) av 

= tf/.(A) + MB)], ( 46 ) 
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where 

A = (a~ + '2ab sin a + 6 2 )*, 

B = (a 2 - 2ab sin a + 6 2 )*, 

(sin (a cos U ) cos [6 sin (a - 0,,)]) = 0. (47) 

Inserting the general formulas equation (46) and equation (47) into 
equation (45), it becomes 

(x m (0;d)x n (D;d))„ = | [MA,) + J Q (BA + J (A 2 ) + J (£ 2 )] 

(x m (0;d)y n (D;d)) av = (48) 

where 

"M = )3[D 2 ± 2Dd(M + 1 - m - n)sma + </ 2 (il/ 4- 1 - m - n) 2 ]* 



.1 



"> = 0[/) 2 ± 2Dr/(m - n) sin a + d 2 (m - w) 2 ]-. 



From equation (48), we can deduce the results 

(x m (0;d)x„(0;d)) llv 

= N{J [(3d(M + 1 - m - n)] + J [/3d(m - n)]\ 

= N[J t) (A„) + J (B )] (49) 

and 

<.c;;(0; d)) = N{JoW(.M + 1 - 2m)] + 1}. (50) 

Then substituting equations (4S) and (49) into equation (44), we com- 
plete the derivation of a normalized correlation function of a signal 
received from a broadside directional antenna array. 
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